Opportunity to Test non-Newtonian Gravity Using Interferometric Sensors with 

Dynamic Gravity Field Generators 
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We present an experimental opportunity for the future to measure possible violations to Newton's 
1/r 2 law in the 0.1 — 10 meter range using Dynamic gravity Field Generators (DFG) and taking 
advantage of the exceptional sensitivity of modern interferometric techniques. The placement of a 
DFG in proximity to one of the interferometer's suspended test masses generates a change in the 
local gravitational field that can be measured at a high signal to noise ratio. The use of multiple 
DFGs in a null experiment configuration allows to test composition independent non-Newtonian 
gravity significantly beyond the present limits. Advanced and third-generation gravitational-wave 
detectors are representing the state-of-the-art in interferometric distance measurement today, there- 
fore we illustrate the method through their sensitivity to emphasize the possible scientific reach. 
Nevertheless, it is expected that due to the technical details of gravitational- wave detectors, DFGs 
shall likely require dedicated custom configured interferometry. However, the sensitivity measure 
we derive is a solid baseline indicating that it is feasible to consider probing orders of magnitude 
into the pristine parameter well beyond the present experimental limits significantly cutting into 
the theoretical parameter space. 

PACS numbers: 



I. INTRODUCTION 



Beyond the original goal of detecting gravitational 
waves, modern interferometric sensors (IFO) (such as 
the future Advanced LIGO (aLIGO) H, the Ad- 
vanced VIRGO (aVIRGO) % the LCGT H, the AEI 
10m [J] interferometer; and the proposed third genera- 
tion detectors such as the Einstein Telescope (ET) @) 
can be viewed as unique pioneering instruments ca- 
pable of measuring induced displacements below the 
1(T 20 - HT 19 m/VHi scale Q. This feature opens 
up myriads of new possibilities for fundamental science. 
In searches for gravitational waves (GW) the signature 
of the induced displacement of the IFO's test mass (TM) 
ranges from the basically unknown to the well-predicted 
but yet undetected. In contrast to this, one may de- 
sign a precision experiment in which the local gravita- 
tional field around the TM of the interferometric sensor 
is modulated at a given frequency on a well-controlled 
manner by a dynamic gravity field generator (DFG). A 
DFG consists essentially of a rotating mass of null odd 
and non-null even moments, from which the quadrupole 
term dominates. In the Newtonian limit, its effect on the 
detector is to add a signal at the even-harmonics of the 
rotation frequency. A well-characterized DFG has the 
potential to provide sub-percent calibration for the GW 
detectors in phase and amplitude as well as to evaluate 
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current theoretical Newtonian noise estimates. This is 
the subject of a separate publication Q. 

It has been shown in the past that devices capable 
of producing gravity field gradients can be employed, in 
conjunction with a suitable displacement sensor, for test- 
ing Newton's 1/r 2 law in the laboratory scale. In the 
first experiment in 1967 Forward and Miller Q used an 
orbiter sensor, originally developed for measuring the lu- 
nar mass distribution, to test Newton's 1/r 2 law in the 
10 cm scale. Weber and Sinsky 0, [Io| used a GW bar 
detector as a sensor, acoustically stressing a volume of 
matter at 1660 Hz, measuring an excess in noise in the 
detector consistent to theory. 

In the 1980s at the University of Tokyo, a series of ex- 
periments were carried through to test violations to the 
inverse square law (ISL) up to a distance of 10 m 

EMI. 

In these studies, the coupling between the dynamic grav- 
ity field generated by a rotating mass, and the quadrupole 
moment of a mechanical oscillator antenna was measured 
as a function of the rotor-antenna separation. With this 
method limits on non-Newtonian gravity were provided 
and a measurement of the Newtonian constant, G, was 
found in agreement with previous experiments [Til [l5j . 

In the 1990s, the gravitational-wave group at the Uni- 
versity of Rome developed and carried out experiments 
[II 0J on the cryogenic GW bar detector, EXPLORER, 
at CERN. A device, rotating at a frequency close to half 
of the antenna's resonant frequency was developed, and 
the resulting dynamic field was measured as a function 
of the source-sensor separation. The results were then 
used to derive upper limits to Yukawa-like gravitational 
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FIG. 1: Current 95%-confidence-level constraints on the 
inverse-square-law-violating Yukawa interactions with A > 
1 cm. (Courtesy of Ref . [la] . Figure(4), with the addition of 
QJ] and [13] data) 



potential violations at laboratory scale. 

The experimental designs cited above consisted in the 
use of a single DFG in conjunction with bar type GW de- 
tectors. In this article we propose a significantly different 
design, exploiting the exceptional bandwidth and sensi- 
tivity of state-of-the-art interferometric sensors. We illus- 
trate the method through the sensitivity of advanced and 
third-generation gravitational wave detectors to empha- 
size the feasible scientific reach. However, it is expected 
that due to the technical details of gravitational wave de- 
tectors, the DFGs based ISL measurements shall require 
dedicated custom configured interferometry. Neverthe- 
less, the sensitivity measure we derive is a solid baseline. 
The concept is based on a null-experiment configuration, 
where a pair of well matched and symmetrical DFGs, ro- 
tating at the same frequency but 90° out of phase, gener- 
ates a null-signal in the absence of violations to Newton's 
1/r 2 law. In the presence of violations, and within exper- 
imental uncertainties, the effect of the two DFGs would 
not cancel, and a measure of such deviations would be 
achievable. 

The most common way of interpreting composition in- 
dependent tests of non-Newtonian gravity (for review see 
e.g. [l8[ and [2l|) is through the Yukawa formalism, 
where an additional so-called Yukawa term (V Y (r)) is 
added to the classical Newtonian gravitational potential 
(V N (r)): 



V(r) = V N (r)+V Y (r) 



-G^ri + ae-r/Aj (L1) 



In this expression G is the gravitational constant, m 
and M are the interacting masses, and r denotes the 
distance between the two point masses. The Yukawa pa- 
rameters of interest are a, denoting the Yukawa interac- 
tion coupling strength, and A, giving the length scale of 
the coupling. For a Yukawa range of A = 0.1 — 10 m, 
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FIG. 2: Schematic of a hypothetical asymmetric two-mass 
DFG configuration. Point masses mi and vri2 are rotating 
with a frequency of fo and radii of r± and r2 from the center of 
rotation, respectively. The DFG center of rotation is placed in 
a distance of d from the mass M, representing the TM center 
of mass. The distance between the DFG masses and the TM 
are h\ and /12, respectively, 6{t) = uot, and only accelerations 
of the TM along the IFO cavity axis are considered. 



the current limit on the coupling strength, a, is in the 
order of 10~ 4 — 10~ 3 for both negative, positive, and ab- 
solute value of a. (see Fig. Q]). The concept designed 
here shows promise to explore deviations from ISL signif- 
icantly below present bounds, maybe even down to the 
order of 10~ 6 if enabled by technology and cost. 

In this work we study the application of future interfer- 
ometric sensors for artificially generated gravity fields by 
a pair of hypothetical DFGs in a null experiment configu- 
ration. We show that a pair of symmetrical and matched 
DFGs cancel each other's effect on the TM at twice 
the rotation frequency in the Newtonian limit. Thus 
we relate the estimated test mass displacement to non- 
Newtonian parametrization of the law of gravitation. We 
also consider the contribution from some of the undesired 
system asymmetries due to production/measurement un- 
certainties. We numerically compute the tolerances of a 
hypothetical experimental setup to measure Yukawa-like 
deviations from the 1 jr 2 law and we estimate a limit on 
the coupling strengths, as, that might be achieved with 
future measurements. 



II. GRAVITY FIELD DYNAMICS FROM A 
TWO-MASS DFG 

Consider a hypothetical DFG, consisting of two point 
masses, mi and mi, separated by a distance of r\ and 
ri respectively from the center of rotation, rotating at a 
frequency of /o = ujq/(2tt). A point mass, M, represent- 
ing the TM, is then placed at distance d (> ri,a) from 
the DFG's center of rotation. 

We first calculate the acceleration of M induced by the 
DFG along the axis connecting M and the center of ro- 
tation of the DFG. This axis corresponds to the optical 
axis of the IFO. Introducing the dimensionless parame- 
ters i?i = ri/d, and Ri = —ri/d, the distance between 
the DFG's i-th mass and the point mass M, hi, can be 
written as 



hi(t) = d Jl + Rf - 2Ri cos 9(t) 



(2.1) 
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where 9(t) — u>ot. The Newtonian potential at the 
position of point mass M is then given by 



where hi = hi (t) is given by Eq. (|2 . 1 1) . In the same way, 
the magnitude of the TM induced acceleration along the 
optical axis is 



(2.2) 



The magnitude of the induced acceleration on M along 
the optical axis can be written as 
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Here Bi(Ri,6) is a geometrical factor 



Bi(R h e) = 



1 — Ri cos ( 



(1 + Rf -2RiCos9) 



3/2 ' 



(2.4) 



As we have shown elsewhere Q , for the case of a much 
smaller lever arm r; than distance d (Ri <C 1), consider- 
ing only the dominant terms in the first few harmonics, 
the induced displacement of a free mass M along the 
axis connecting it to the DFG's center of rotation can be 
approximated as 
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where AA n describes the n-th multipole moment of 
the DFG's mass distribution. We observe that using 
Eq. (|2.5p and in the case of a DFG mass distribution 
invariant under rotation by pi, all odd moments van- 
ish and the induced displacement is dominated by the 
quadrupole moment M2 at twice the rotation frequency. 
We note, that the solution represented in Eq. (I2.5[) is valid 
for a free body, which, in case of a suspended TM is a 
good approximation for frequencies well above the eigen- 
frequencies of the suspension (typically around 0.1 — 1 Hz 

In the same way as above, we calculate the field dy- 
namics arising from a Yukawa-like potential perturbation 
to the classical Newtonian field. In view of the fact that 
the potentials are additive and all operations to calculate 
the acceleration are linear, the acceleration of point mass 
M due to the two potential terms will also be additive. 
Using Eq. (jl.ip . the Yukawa perturbation to the Newto- 
nian potential at M due to the DFG can be expressed 
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where 9 = 9{t) = ui n t, and fi(Ri,9,X) is a function of 
the length scale of the Yukawa coupling, A: 
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Analytical integration of Eq. (12.31) and Eq. (|2.7[) cannot 
be easily formulated without using some additional level 
of approximation. For this reason, we treat the prob- 
lem numerically, facilitating the treatment of some of the 
experimental uncertainties. We have modeled both the 
Newtonian and the non-Newtonian dynamics on the TM 
while taking into account many of the DFG fabrication 
and measurement procedural uncertainties. 

In this work the TM is approximated as a driven, 
damped pendulum. Let the pendulum's longitudinal 
eigenfrequency be uj p with quality factor Q. Casting the 
differential equations of motion (Eq. (|2.3p and Eq. (|2.7p ) 
in the Laplace frequency domain, the transfer function 
relating the acceleration a N ' Y to the induced displace- 



ment cc N ' Y is 



„n.y 
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(uj p /Q)s + u 2 



(2.9) 



where s = iu (the upper index N, Y denotes that the 
same expression applies to both Newtonian and Yukawa 
case). 

We used Monte Carlo simulations to model a null- 
experiment design using DFGs to measure hypothetical 
violations to the classical inverse square law. To do this, 
we first numerically computed the induced acceleration 
due to the Newtonian and non-Newtonian terms using 
Eq. (123)1 and Eq. (|2~7|) . These results, generated in the 
time domain, are then mapped in the Fourier-domain to 
generate an acceleration spectrum, which is then filtered 
through the transfer function presented in Eq. (|2.9p . The 
displacement contributions at the different harmonics, for 
different sets of parameters and uncertainties are then 
computed and analyzed. 

To assess the short and long term feasibility of the pro- 
posed method, the simulated results are compared to the 
displacement sensitivity of future most sensitive interfer- 
ometric GW detectors. Fig. [3J shows the design sensi- 
tivities for aLIGO, aVIRGO, the ET (for a reference to 
these, see Q), LCGT 0], and the AEI 10m detector j||. 
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FIG. 3: The design displacement sensitivity of aLIGO (thin 
black), aVIRGO (dashed), the ET (dot-dashed), LCGT (dot- 
ted), and AEI 10m (thick green). 




FIG. 4: Null experiment geometry for the measurement of a 
possible deviation from Newton's 1/r 2 law. Two conceptual- 
ized DFGs are placed in-line with one of the IFO's arms at a 
distance of d 1 and d 11 from the TM. The DFGs positions and 
quadrupole moments are chosen such that when operated at 
a relative phase of /3 — it/2 the net displacement at 2/o will 
largely depend on the non-Newtonian coupling strength a. 



III. HYPOTHETICAL INVERSE SQUARE LAW 
VIOLATION MEASUREMENTS WITH DFGS 

In order to mitigate the effect of the IFO's calibration 
uncertainty, two DFGs rotating at the same frequency 
(/o), but out of phase (f3 — 7r/2), are set up to generate 
a null effect on the TM in the Newtonian limit at twice 
the rotation frequency. For this reason, the uncertainties 
of the proposed measurement, whose setup is shown in 
Fig- HI, will solely depend on the production and proce- 
dural uncertainties of DFGs. 

For example, similarly to the one proposed in [7| as 
a detector calibration tool, a single hypothetical DFG 
may consist of an Aircraft Grade (6Al/6V/2Sn) Tita- 



nium disc 60 cm in diameter and 10 cm in height. The 
disc has two cylindrical slots, 40 cm apart (r\ 2 = 20 cm), 
which can hold different materials. Tungsten cylinders, 
16.1 cm in diameter, can serve as rotating masses, with 
rn\ 2 = (pw ~~ PTi)V = 30 kg effective mass each, where 
pw and pTi are mass densities of tungsten and titanium 
respectively, and V is the volume of a cylinder. 

In the ideal case of symmetrical and matched DFGs, 
no effect on the TM is observed at 2/o in the absence of 
1/r 2 violations. In the classical sense, this condition is 
satisfied if the DFG radii and placement follow the 

r l i,2 = Vvr{, 2 (3.1) 
d 11 = ^/rjd 1 

scaling laws, where indices I and II are associated with 
the two DFGs, and r\ = m^/m* 2 is the mass ratio (see 
Fig. 01). Note, that if we take into account the cylin- 
drical geometry of the TM, the scaling laws presented in 
Eq. (I3.1[) will differ by O(0.1%), where the corresponding 
r l i 2 and d 11 will still be analytically calculable. 

In non-Newtonian dynamics, using the Yukawa formal- 
ism, limits on the parameter alpha are given as a function 
of the Yukawa range, lambda. First we address an ideal 
case, where the parameters of both DFGs are known ex- 
actly. As an example, let's assume the first DFG, with 
parameters described above, positioned at d 1 — 2 m from 
the center of mass of the IFO TM. Specifications and po- 
sition for the second DFG are determined by Eq. (|3.ip . 
On Fig. [5] we plot the interaction range for which the 
TM displacement at 2/o is maximal, A max , as a function 
of i] (continuous line). On the same plot, the dashed 
line shows the corresponding TM RMS displacement at 
2/o = 10 Hz due to the Yukawa term. 

Using the above case and 77 = 2, Fig. [5] gives an inter- 
action range of A max = 0.64 m, corresponding to a max- 
imum displacement of the TM due to the Yukawa term 
with an RMS value of x Tms = (6.8 x 10~ 16 m) x \a\ . 
For the case of the AEI 10m detector with a noise floor 
of n = 1.5 x 10~ 17 m/VHz at 10 Hz, and an integra- 
tion time of T = 1 day, a limit of \a\ ~ 2.3 x 10 -4 can 
be provided with a signal to noise ratio (defined as the 
ratio of the RMS signal to the displacement noise spec- 
trum density integrated for a time T) of SNR = 3, given 
the basic physical limitations described later in this sec- 
tion. In general terms, for an arbitrary noise floor n, and 
integration time T, the signal to noise ratio scales as 

^^-( '■ 5 " c ; m/ i H"' m 

/ Xrms \ ( |oj \ 
I 6.8 x 10" 16 m M 2.3 x 10" 4 J 
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FIG. 5: Simulation results for two DFGs positioned in a 
null experiment. The continuous line represents the plot of 
A max (?7) expressed in d 1 units, where A max denotes the Yukawa 
interaction range, for which the 2/o RMS displacement due 
to the Yukawa term is maximal. The corresponding RMS dis- 
placement of the TM (using d 1 = 2 m, r 1 = 20 cm, m 1 = 30 kg, 
and fo=5 Hz), expressed in \a\ units, is shown by the dashed 
curve. 



A. Measurement Uncertainties 

There arc two fundamental classes of effects that 
limit our capability of measuring an inverse-square-law- 
violating gravitational force within the experimental con- 
cept. The first one is due to the noise level of the IFO, 
which effect can be suppressed by using sufficiently long 
integration times, as discussed in Sec. IIII1 In this sec- 
tion we discuss some of the effects from the second class, 
which are due to the finite measuring and manufacturing 
precision when setting up the null-measurement. 

Errors due to conceivable machining precision and pro- 
cedure, uncertainties related to measurement of param- 
eters of the final DFG products as well as their final 
placement relative to the TM and to each other will all 
infer an uncertainty in the induced displacement mea- 
surement, and propagate an error into the measurement 
of the non-Newtonian parameter a. In case of purely 
symmetric DFGs positioned in a null experiment config- 
uration, and taking into account only the Newtonian field 
component, the induced displacement of the TM will be 
dominated by a 4/o term. In the presence of asymme- 
tries, the displacement will also have components at the 
odd harmonics and the classical terms at 2/o will not 
cancel each other completely. This means that the New- 
tonian peak due to this imperfect cancellation will appear 
at 2/o, and will contaminate a potential Yukawa-signal 
from the very beginning of the measurement. 

First, we calculate the achievable lower limits on the 
Yukawa strength parameter, \a\, in case of ideal DFGs 
with no parameter errors. We chose the parameters of the 
first DFG to be m l = 30 kg, d l = 2 m, and r l = 0.2 m, 
and maximized the integration time of the measurement 



at T = 10 7 s (~ 4 months). Because of the slight depen- 
dence of the TM response on r\ (see Fig. \5\), and taking 
into account manufactural limitations, we set r\ to 2. We 
chose a conservative limit of SNR — 8 as the condition 
for detection. The optimal DFG frequencies that maxi- 
mize the SNR in terms of |a| for the different detectors 
(see Fig. along with the lowest achievable limits on 
the strength parameter, \a*\, are given in Table HI Re- 
sults on \a\ as a function of the Yukawa scale parameter, 
A, for five different IFOs are provided in Fig. [5] Note, 
that these results are also valid for the case, when the 
measurement setup errors are kept low enough to allow 
the detector noise levels to be the main limiting factors. 

TABLE I: Results of the numerical calculations, excluding 
DFG-related experimental uncertainties, for achieving abso- 
lute limits on \a\ in case of maximal detector response due 
to the Yukawa term at 2/o for an integration time of 

T = 10 7 s ~ 4 months. The DFG parameters used for gener- 
ating the results are m 1 = 30 kg, S = 2 m and r 1 = 0.2 m. 
Amax corresponding to \a* \ is A max = 0.64 m in all five cases. 
The results on optimal DFG frequencies for ET and AEI 10m 
are limited by the uncertainties of the IFO sensitivity models 
at low (< 10 Hz) frequencies. 



T = 10 7 s 


aLIGO 


aVIRGO 


ET 


LCGT 


AEI 10m 


fo 


11 Hz 


20 Hz 


<5Hz 


26 Hz 


<5Hz 


Kl 


1.2 x lCT 6 


1.6 x HT 6 


1.9 X 1CT 7 


3.4 x 1CT 6 


5.6 x 10~ 5 



In order to study the feasibility and real-life limita- 
tions of null-experiment concept, we numerically studied 
the effect of uncertainties associated with the system pa- 
rameters included in the simulations (rotation frequency 
(fo), mass (m) and radius (r) of each DFG, as well as 
the DFG positions relative to the TM and to each other 
(d and /?)), thereby mimicking a realistic construction 
and measurement procedure. To do this, a large number 
(N) of hypothetical setups were Monte Carlo simulated 
with the DFG parameters normally distributed around 
a mean. The mean of the parameter distributions (e.g. 
{m}) — 30 kg, (d l ) = 2 m and (r l ) = 0.2 m) are chosen 
such as to maximize the response of the IFO in terms 
of | a |, while taking into account the technically achiev- 
able range of values for each parameter, as well as safety 
issues. 

We again set ij to 2 for all setups, similarly to the ideal 
case. Following Eq. ()3.1[) . r\ in this case determines the 
means of the distributions corresponding to parameters 
of the second DFG. The means of DFG operational fre- 
quencies were chosen to be equal to the values given in 
Table |TJ The low optimal DFG frequencies allow the ef- 
fective masses of DFG fillings for the first DFG (rotating 
closer to the IFOs' vacuum chamber) to m 1 — 30 kg. 
The highest corresponding mechanical energy of a DFG 
(to 11 = 60 kg; r n =V2x 0.2 m; f = 26 Hz) would still 
be less than 0.02% of the energy proposed in a similar ex- 
periment described in (23|. Note, that as an alternative 
to the DFG geometry we proposed in Sec. IHIja rotating 
mass design similar to the one suggested by [23| can also 
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be used in the measurement. 

Naturally, we are seeking to achieve the highest pre- 
cision measurement currently possible or feasible in the 
near future, and we let the parameter uncertainties vary 
to values corresponding to state-of-the-art limits with 
the investment of plausible amount of work and finances. 
Our goal is to find, using Monte Carlo simulations, the 
limits on \a\, allowed by using a two-DFG configuration 
with parameter uncertainties. 

Each parameter associated to hypothetical setups were 
randomized independently, from a Gaussian distribution, 
using the same variances for corresponding parameters 
of different setups. We used a chosen set of variances for 
the parameter distributions (discussed in Sec. IIII Bp that 
are within the limits of current state-of-the-art machin- 
ing and measurement technologies, and at the same time, 
minimize the Newtonian peak we get at 2/o due to im- 
perfect cancellation. We generated N — 1000 two-DFG 
setups, and determined the 95% quantile of the induced 
Newtonian displacement values at 2/o (from now on, de- 
noted by Xg 5 ). Note that this quantile x$ 5 is completely 
independent from the detector noise and thus from the 
measurement integration time. We then, for each A, cal- 
culated the \a\ value for which the resulting Yukawa peak 
at 2/ equals Xg 5 . The limits on \a\ defined this way, ver- 
sus A, are shown in Fig. [6] 

Using a two-DFG setup with allowed parameter uncer- 
tainties given in Sec. IIII B[ the limiting factor in case of 
the LCGT and AEI 10m IFO is its relatively high noise 
level at 2/o. However this means, that we can toler- 
ate higher parameter uncertainties compared to the ones 
given in Sec. IIII B[ when using LCGT and AEI 10m and 
an integration time of 4 months. In case of the other 
three IFOs (aLIGO, aVIRGO, and ET) the limiting fac- 
tor is the imperfect cancellation of the Newtonian peaks 
at 2/ due to parameter errors, when considering a 4 
months long integration time. Therefore in order to reach 
the best possible \a\ limits for the chosen set of parame- 
ter uncertainties, we can use a reduced integration time 
of O(10 4 — 10 6 s), depending on which IFO we choose. 
The reach of the allowed exclusion region suggests that 
the displacement sensitivity of any dedicated IFO for the 
null-measurement should have comparable displacement 
sensitivity to advanced interferometric gravitational wave 
detectors. Custom configured interferometry could pos- 
sibly have several advantages such as double DFG com- 
patible suspensions, being able to reduce the distance be- 
tween DFGs and the IFO test mass, and thus to obtain 
better measurement results. 



B. Measurement feasibility 

Monte Carlo test results presented in Sec. IIII Al were 
provided by using a chosen set of DFG parameter uncer- 
tainties, that gives the best coverage of Yukawa parame- 
ter predictions given by the " t" model (see [2l[ for more 
details), and that are still within the feasible measuring 
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FIG. 6: Predictions for the allowed lower limits on Yukawa pa- 
rameter | a | via laser interferometric method proposed here. 
The grey shadings show the parameter regions excluded by 
previous measurements [Isl . [24[ that were sensitive to neg- 
ative, positive, and absolute alpha values only. The thick 
black line corresponds to the limit achievable with an ongo- 
ing experiment, denoted by UW-UCI (24|. The red curves 
show the results for various IFOs in case of infinite DFG ma- 
chining precision, using an integration time of T = 10 7 sec. 
The IFO+PA curve gives the 95% confidence limits for all 
IFOs but LCGT and AEI 10 m, wh en allowing DFG parame- 
ter uncertainties given in Sec. IIII Bl The limiting factor in the 
measurement is the finite integration time when using LCGT 
and AEI 10m, and the finite machining precision in case of all 
the other IFOs. The coloured background areas correspond 
to different theoretical models presented in [2ll |. that predict 
Yukawa parameters within the areas. Our tolerable parame- 
ter uncertainties were chosen such as to give the best coverage 
of the parallelogram corresponding to model "t", while taking 
into account the feasible parameter adjustment precision. 



and machining precision. An optimization of tolerable 
DFG parameter uncertainties in terms of measuring an 
ISL-violating gravitational force can be given with the 
Fisher matrix method [IBj], however this analysis is be- 
yond the scope of this proof-of-concept paper. In this 
section we discuss the DFG parameter uncertainties that 
correspond to the IFO+PA curve given in Fig. [6l Note, 
that we assumed that the parameter uncertainties are 
kept at the same levels throughout the whole integration 
time of the Yukawa measurement. We also discuss the 
possible methods in measuring and machining that allow 
us to reach the chosen DFG parameter precision. 

In the conceptual experimental setup considered the 
two DFGs are rotating at fo = 11 Hz, 20 Hz, 5 Hz, 
26 Hz, and 5 Hz using aLIGO, aVIRGO, ET, LCGT, 
and AEI 10m, respectively. In all cases the two-DFGs 
in a setup have to be in a relative phase of /3 = n/2. 
In our simulations we chose the following measurement 
uncertainties: 

• Uncertainty in the frequency of the first DFG (ab- 
solute measurement): Sfi = 10~ 6 Hz 
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• Uncertainty in the frequency of the second DFG, 
relative to the first DFG (relative measurement): 
8f 12 = 1CT 9 Hz 

• Uncertainty in the initial phase difference be- 
tween the two DFGs (absolute measurement): 
Sp = fCT 8 rad 

Frequencies and the relative phase between the differ- 
ent DFGs can be finetuned to their desired value by cal- 
ibrating and locking the DFGs' rotation signal with an 
accurate and low-phase- noise clock (e.g., good Cesium 
or GPS), also used for time-stamping the IFO data. A 
sophisticated DFG frequency (and phase) control system 
can be constructed allowing nanosecond timing precision. 

In our simulations we used tungsten filling masses of 
30 kg for the first, and 77 x 30 kg for the second DFG. 
The masses had the following permissible uncertainties: 

• Uncertainty in the exact mass of one of the tung- 
sten cylinders in the first DFG (absolute measure- 
ment): Smi = 5 x I0~ 5 kg 

• Uncertainty in the mass of any other DFG cylinders 
of the two-DFG setup (relative measurement). The 
mass of one of the tungsten cylinders built into the 
second DFG has to be 77 times the mass of the one 
in the first DFG, where we proposed to use r\ — 
2. This way we can use a relative measurement, 
and thus, we have to adjust the mass of one of 
the cylinders in the second DFG to the sum of the 
masses in the first DFG. We adjust the mass of the 
second tungsten cylinder in the second DFG to the 
mass of the first cylinder in the same DFG. The 
mass of the second cylinder in the second DFG can 
be then finetuned to the mass of the first cylinder 
within the same DFG. The tolerable uncertainty of 
these relative measurements and mass adjustments: 
8m 12 = 1 0~ 6 kg 

For the DFG mass measurements, one can use scales 
such as the Sartorius ME415S scale (26|, that can deter- 
mine the absolute value of the masses with an uncertainty 
of ~ I0~ 7 kg with 410 g weighing capacity. Using state- 
of-the-art mass comparators allows us to reduce the rel- 
ative measurement (mass-to-mass comparison) error to 
~ 2 x I0~ 9 kilograms (e.g 27]). A complete tungsten 
cylinder can be built from an assembly of several smaller 
disks, allowing us to reduce uncertainties with precision 
scales having lower weighing capacity. These smaller and 
lighter components of the tungsten fillings are also easier 
to be manufactured and handled. A combination of pre- 
cision absolute scales, mass comparators, and fine mass 
adjustment through abrasion leads to precise mass stan- 
dards and also allows manufacturing of matched mass 
pairs beyond the load limit of the absolute scale. 

Arm lengths and DFG positions have to be set to def- 
inite values (in our example case r\ 2 = 0.2 m, r{ J 2 = 
\/2r\ 2 , d 1 = 2 m, and d 11 — ^/2d l ) to achieve precise 



cancellation of the Newtonian component. The corre- 
sponding parameters have the following uncertainties: 

• Uncertainty in the distance between the first DFG's 
rotation center and the TM's center of gravity {d l 
absolute measurement): 8d\ = I0~ 2 m 

• Uncertainty in the distance between the two 
DFGs relative to d l (relative measurement): 
6di2 = 10- J m 

• Uncertainty in the length of one arm in the first 
DFG (r 1 absolute measurement): 8r\ = I0~ 3 m 

• Uncertainty of in the equality of arms within one 
DFG (fru^)) an d the uncertainty in the length 
of one arm of the second DFG relative to r 1 (8r 12) 
(relative measurements): 8r\\^2 = <5^i2 = 3 x 
I0~ 8 m 

Distance-like quantities (r-s and d-s) can be measured 
using interferometric methods with accuracy significantly 
surpassing the wavelength of the laser light. The exact 
technique of the measurement greatly depends on the 
DFG placement, accelerating and controlling configura- 
tion, and is beyond the scope of this article. 

Setting the arm lengths within one DFG to be equal 
with a precision of 3 x 10~ 8 m is challenging, but possibly 
feasible. The uncertainty of equivalence of arms within 
one DFG will be in the order of the equivalence of masses 
within the DFG, and the uncertainty of Aii being zero. 
The feasible precision of mass-to-mass comparison can 
be as low as ~ 2 x I0~ 9 kg [27}. In case of a non-zero 
A4i, the axis of rotation of a spinning DFG will subject 
a periodic driving force (that is proportional to ~ /q), 
causing the support of the axis to vibrate. The ampli- 
tude of the vibration can be measured with high-precision 
interferometry, and the placement of the center of rota- 
tion can be adjusted until the vibration amplitude is zero 
within the measurement precision. The adjustment can 
be carried out using much higher DFG frequencies than 
the values given in Table HI allowing an amplification of 
DFG vibrations in case of a non-zero AA\. 

The non-Newtonian and the classical gravitational po- 
tentials scale differently with d, the distance between the 
TM and the center of rotation of a DFG. Using the two 
DFGs positioned in a null experiment configuration, and 
measuring the TM displacement spectrum in case of dif- 
ferent d-s, can also help us to determine measurement 
and setup inaccuracies other than errors of d. 

It is possible that periodically varying one or more of 
the tunable parameters (e.g., distance or phase) will al- 
low precise, in — situ tuning of the cancellation neces- 
sary for the null- measurement. Also, an additional effect 
can be exploited to estimate the dominating error in the 
two-DFG setup parameters. By applying the parameter 
uncertainty tolerances given in this section, the TM's in- 
duced displacement at /o due to imperfect cancellation of 
the Newtonian term will be orders of magnitudes above 
the detector noise level and the Yukawa term at /q, and 
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thus will be detectable. By monitoring this signal at fo 
we can give a rough estimation on the dominating error 
in the parameters of the two-DFG setup. 

We investigated the projected contributions of the dif- 
ferent parameter uncertainties, relative to each other, to 
the uncertainty of an \a\ measurement. To do this, we 
carried out a series of Monte Carlo simulations, where 
in each simulation, we set the relative uncertainty of a 
chosen parameter to 10~ 6 , and kept the uncertainties 
of all the other parameters at zero level. We simulated 
N = 1000 two-DFG setups, and calculated the relative 
error of \a\ obtained from the N = 1000 hypothetical 
measurements of |a|. In order to characterize the rela- 
tive dominance of the different parameters compared to 
each other, we normalized the measured relative errors 
of | a | with the highest such relative error value we ob- 
tained. This way, our final results were basically indepen- 
dent from the 10 -6 relative uncertainty chosen for each 
individual setup parameters. The resulting relative dom- 
inances of the different parameter uncertainties in the |a| 
measurement are given in Table \H\ Our results show that 
the dominant parameter - for which the allowed uncer- 
tainty, 5di2, should be kept as low as possible - is the 
distance between the two DFGs, relative to d 1 . 

TABLE II: The relative dominance (RD) of parameter uncer- 
tainties in terms of their contribution to the uncertainty of 
a | in a hypothetical measurement. We carried out TV = 1000 
Monte Carlo trials for each uncertainties, setting the chosen 
uncertainty to a 10 -6 relative value, and while setting all 
the other uncertainties to zero. The relative dominance for 
each parameter uncertainty was obtained by calculating the 
relative uncertainty of \a\ measured from the N trials, and 
normalizing it by the largest such relative uncertainty value 
(thus we get a relative dominance of 1 for the most dominant 
parameter uncertainty). 

Uncertainty <5di2 5d\ &fi2 <5>"i2 <5j"ii,22 
RD 1.0 0.52 0.37 0.32 0.28 

Uncertainty Sfi 5r\ 5m\2 dmi 5j3 
RD 0.26 0.26 0.217 0.127 2.3 x 10 -6 



In this subsection we addressed the fundamental 
sources of measurement errors and the uncertainty limits 
for each parameter of interest. We also pointed out, that 
uncertainty values for DFG parameters used in our pro- 
posed experimental setup are within the limits of current 
state-of-the-art machining and measurement technolo- 
gies. For a more realistic application a full dynamic finite 
element simulation of the mass distributions must be per- 
formed. The expansion and stress factors of the DFGs 
under prolonged operating conditions must be modeled 
and simulated, then subsequently measured and taken 
into account. Practical limitations due to the complex 
mechanical and detector dynamics/geometry around the 
test mass and other second order error sources need to 
be investigated, via experiments and more sophisticated 
simulations as they might prevent us from reaching the 



full potential of DFG device based measurements. For 
example, in practice material inhomogeneities can give 
rise to a substantial uncertainty that could be minimized 
by using NDT-rated materials [28[ and also reduced by 
different geometries. Alternative, more simplified DFG 
geometries different from the one proposed here (e.g. ro- 
tating rods) could also possibly lead to higher measuring 
and machining precision, and therefore should be stud- 
ied carefully. All these issues will be the scopes of future 
studies. 



IV. SAFETY 

Significant kinetic energy (i.e. ~ 100 kJ) is stored in 
the DFG once it rotates, therefore crucial safety consid- 
erations must be addressed. There are two major points 
of failure management to be concerned with. 

(a.) The vacuum chamber of the DFGs must be made 
strong enough to withstand the damage of an accidentally 
disintegrating disk. This is the standard solution for high 
speed gyroscopes. 

(b.) For added security, the gap between the inner wall 
of the vacuum chamber and the outer edge of the rotat- 
ing disk could be kept relatively small. In the event of 
an incident where the DFG's material starts to yield or 
its angular acceleration is uncontrolled the disk will ex- 
pand radially touching the sidewall and may slowly stop, 
potentially preempting some of the catastrophic failure 
modes. Adding a heavy ring in contact with the inner lat- 
eral wall of the vacuum chamber (similarly as described 
in would absorb a huge angular momentum of the 
DFG and dissipate its rotation energy via the friction 
against the lateral wall. 

These conditions can be met using Finite Element 
Analysis (FEA) aided design, in-house destructive test- 
ing of sacrificial parts and relying only on the best base 
materials. 



V. CONCLUSION 

Two DFGs in a null experiment setup in conjunction 
with an interferometric sensor allow for studies of com- 
position independent non-Newtonian contribution to the 
classical gravitational field in the meter scale. Simula- 
tion results presented in this paper indicate that by tak- 
ing advantage of two matched DFGs and the relatively 
large bandwidth and high sensitivity of state-of-the-art 
interferometric sensors there is an exciting opportunity 
to explore a below the current limit in the A ~ 0.1 — 10 
meter Yukawa range addressing standing theoretical pre- 
dictions. 

For the proof of concept study, we chose a conserva- 
tive 2 meter distance for our device from the test mass 
of the interferometric detector and used the sensitivity of 
advanced and third-generation gravitational wave detec- 
tors as the realistic baseline. We intend to note here, that 



9 



putting the devices considerably closer to the test mass 
can yield orders of magnitude better limits in ISL vio- 
lation parameters through the characteristic length scale 
changes for the Yukawa case. 

There are many practical details that still need close 
attention when designing and manufacturing a practical 
device. Finite element analysis of the DFGs and subse- 
quent experimental studies are necessary to completely 
understand the stresses the DFG is subjected to. The 
DFGs need to be enclosed in separate vacuum chambers 
seated on seismic attenuation stage(s). Prototype design 
and test will be necessary to balance the disk and test 
vibration control. Direct gravitational coupling into the 
test mass suspension and seismic isolation system must 
be carefully designed, simulated and mitigated. Any 
structure supporting the test mass must be placed far 
away (i.e., tens of meters) from the double DFG assembly. 
Other mostly practical problems, such as safety, can also 
be solved as was shown in past applications/experiments 
that have used rapidly rotating instruments (for exam- 
ples see references in Sec. U). It is also likely that ded- 
icated and custom designed interferometric sensors will 
be necessary and the direct use of advanced interfero- 
metric gravitational wave detectors will not be ideal for 
technical reasons. 

In spite of the promising estimates one should not un- 



derestimate the technical difficulties imposed by (1) the 
required design precision, (2) the manufacturing of the 
DFGs and (3) the unknown, but experimentally possible 
to study, complex nature of the coupling routes of the 
dynamic gravity fields of the DFGs into the complicated 
mechanical structure of the IFO test masses. 
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